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Motivation 


Introduction and Motivation 
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QCD phase diagram 


Schematic figure: 


Temperature T 


` Quarkyonic 
Matter 


Baryon Chemical Potential us 
(Fukushima, Hatsuda, Rept. Prog. Phys. 74 (2011) 014001) 


e Finite temperatures --- Lattice simulation works well ! 
e Finite density --- No ab initio calculations. 
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Lefschetz-thimble path integral 


Where is the finite-density QCD? 


Neutron star 

e Cold and dense nuclear matters 

@ 2m, neutron star (2010) 

e Gravitational-wave observations (2017~) 
Heavy-ion collision 


e Low-energy scan of heavy-ion collisions are planned and run in 
many facilities (RHIC, SIS, J-PARC) 


Reliable theoretical approach to equation of state must be developed! 


Z(T, u) = foa Det (P(A) + m) exp —Sym(A). 
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Silver Blaze problem in finite-density QCD 
QCD & |QCD| 


ZQCD — [PA (dety,D,) gr P. ZiQCD| = [ra |dety,D,| B. 


At u = 0, these two are the same! But, 


ng 
QCD 


|QcD| 


m /2 my/3 "H 


|QCD] experiences the phase trans. at u = m4/2 ~ T0MeV. 
In QCD, the state = QCD vacuum for u < my/3 ~ 300MeV. 
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Motivation 
Path integral with complex weights appear in many important 
physics: 


e Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


e Gauge theories with topological 0 terms 


9 Real-time quantum mechanics 


Oscillatory nature hides many important properties of Z. 


Purpose of this talk 


Let us grasp the significance of the silence of Zocp in the region 
lia > ed 2- 
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Lefschetz-thimble integrations 


Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) — EL (= + 2 


Complexify the integration variable: z = x + iy. 


Im(¢) 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
d 3 
Ai(a) — > No [ Pi (5 +a) : 


ng: intersection number of the steepest ascent contour Ke and R. 
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Figure: Lefschetz thimbles J and duals K (a = 1e®", —1) 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent" cycles Js: (classical eom S'(z,) = 0) 


d^z e-$60 = n. | d^z e- 5), 
i 2j, 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


d (20) tim. ai) = ca} . i H (8 


No: intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = zo 3). 


[Witten, arXiv:1001.2933, 1009.6032] 
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Tips: Airy integral & Airy equation 


Let us consider the "equation of motion": 


dz d ei ra) 


= 0. 
2n da 


This is nothing but the Airy equation: 


Two possible integration contours 
— Two linearly independent solutions of eom. 
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Sign problem & One-site Hubbard model 


Simplest model study to understand the sign problem 
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Sign problem & One-site Hubbard model 


One-site Fermi Hubbard model 


Hubbard model (&: annihilation op., fj = ĉtê): 


= -t o» Creag UM ifia -uY fn; + fry i). 


(1,9) o 


This is too difficult as a toy model. 
Turn off the hopping, t = 0: 


A = Uff, — u(fu + f). 
Fock state gives the partition function immediately! 


Z=tr [exp -B| — 14 2e8^ + eBQu-U) 
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One-site Fermi Hubbard model 


The number density is given as 


18 Deh + aP a UP) 
B m In Z ui 1 + Jebu -+ e82u—-U) ` 


In the zero-temperature limit, 


Question 


How can we describe these step functions using the path-integral 
formalism? 
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Sign problem & One-site Hubbard model 


Path integral for one-site model 


The partition function reads 


N p? 
avi f Teo (32585) 
pr dvd, exp — (>: Viet (vest — ef Mu) : 


with the antiperiodic boundary condition iw., = — 41. 
ọ is an auxiliary field for the number density: 


(A) = (yt, em Get) = — > (pp). 
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Sign problem & One-site Hubbard model 


Integrate out fermions 


Fermion spectrum: 


N 
Alp, ui) = 1 — eO DN exp E » ig, + u + U/2). 


It couples only to the Matsubara zero mode of y. 


The path integral reduces to the zero-mode integral: 


J B i U/2)\2 .— 8? /2U 


fermion det. 


Integrand has complex phases causing the sign problem. 
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Sign problem & One-site Hubbard model 


Classical action & flow equation 


Effective action: 


S(z)- faan (reos (iie 7). 


It has a charge conjugation symmetry: 


This ensures the manifest reality of observables. 
(YT, Nishimura, Kashiwa, PRD 91, 101701(R) (2015)) 


The flow equation reads 


de £.. 218 exp B (—iz + u + £) 


dt U^ 1+ exp @ (-iz+p+ 9) 
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Sign problem & One-site Hubbard model 


Flows at u/U < —0.5 and u/U > 1/5 


Re(z) 


Figure: Flow at u/U 


2 


| dz e 8), 
Se 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 
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Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Flows at —0.5 < u/U < 1.5 


025 05 0.75 
Re(z) 


Figure: Flow at u/U = 0 


-0.5 


-0.75 


Behaviors of the flow become totally different! 


u/U + 1/2. 


~N 


Furthermore, saddle points lie around Im(zm)/U 


This value is far away from n = Im (z)/U=0, 1, or 2. 
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Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Curious incident of n in the one-site model 


We have a big difference bet. the exact result and naive expectation: 


05 0 05 1 15 — uU 


CL reproduces the naive one! (vr, Hidaka, Hayata, arXiv:1509.07146, 1511.02437) 


This is similar to what happens for QCD and |QCD|. 
p/U = —0.5 & p=m,/2. 
(cf. Monte Carlo with one-thimble approx. gives the naively expected results. 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 
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Semiclassical analysis 


Semiclassical study of the sign problem 
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Semiclassical analysis 


Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


u[U x 1/2 


U 
Ea cR By + 2xml. 


At these solutions, the classical actions become 


Sox -£ (B +3)’, 
Re (Sin — So) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 
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Semiclassical analysis 


Semiclassical partition function 
Using complex classical solutions Zm, let us calculate 


oo 
Za t= , p nm. 


Tn—-—oo 


This expression is valid for —1/2 < u/U < 3/2. 
This is calculable using the elliptic theta function: 


Za 


K 


= 1 dcc 
-S [149 2 H+ \ q-2n?m?/8U 
e + 2, cos 27m (5 + 2 e 


Lo 90 Lad 1 —2n? [BU 
e ^ (n(5*5).« ) 
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Semiclassical analysis 


Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e o% 
. è 
r . Re 
° e 


Semiclassical study gives the correct transition! 


" 2 (1<p/U « 3/2), 
Nel = 39, Za 1 (0 « p/U <1), 
Pap 0 (-1/2<p/U « 0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Semiclassical analysis 


Important interference among multiple thimbles 


Let us consider a "phase-quenched" multi-thimble approximation: 


Za = $e | = 6750905 (0, e 0797). 


e Lee-Yang zeros cannot appear at u/U — 0, 1. 
e One-thimble, or "phase-quenched" , result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Ünsal, PRL 115 041601) 
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Semiclassical analysis 


Numerical results 
Let us remind the result for n: 


— 1-thimble 
---3-thimble 


We need five thimbles to describe the transition at GU = 30. 


Question 
How many thimbles are necessary for practical purpose? | 
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Semiclassical analysis 


Necessary number of complex classical solutions 
Criterion to neglect the Lefschetz thimbles Jm for large |m]: 


m m-—1 
A T 
gm 


cl 


e 


Assuming |m| >> 1, we can write this as 


2| exp — Sml <e 
Za ii 


Here € < 1 is a controlling parameter of an error. 
We can solve this inequality with respect to m: 


1 
er (en) 
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Semiclassical analysis 


Necessary number of complex classical solutions 


If the sign problem is sever, 05 (geev) ~ 4/2 e-8U/8. 
The criterion gives 


2m 


If the sign problem is mild, 42 ~ 05 (m, eth pur EE. 


We only need 


Necessary number of thimbles depends on the region which one 
wants to study. 

For most severe cases, it is proportional to GU /2z(— oc). 

(YT, Hidaka, Hayata, 1509.07146) 
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Semiclassical analysis 


Numerical results 
Results for GU = 100: 


n 

.5 
gu =100 P| Pod 

exact pp LO D 

--- 7-thimble | bol | - ` 
—17-thimble| L| | k 

I 1 \ 

1 1 

Ws og 

1 1 

i} 1 

1 1 

t i 

d V 

4 \ 


Necessary numbers of Lefschetz thimbles ox BU /2r. 
(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Complex Langevin and Lefschetz thimbles 


Relation between complex Langevin and Lefschetz thimbles 
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Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz,(0) OS 


ag 7 gz 00) + V(O). 


0: Stochastic time, 7: Random force satisfying 
(n(8)n(0')), = 20(0 — 0"). 
Ito calculus shows that 


5 (Olen(0)) = A(O" (z,(8))), — (O'(Zn(9))S"(2n(9)))n- 


If the I.h.s becomes zero as 0 — oc, this is nothing but the 
Dyson-Schwinger eq. 
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Complex Langevin and Lefschetz thimbles 


Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(05) 7 7 374. | ase oq), 


c o 


in d; € C. To reproduce physics, d, = (Ko, R). 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Complex Langevin and Lefschetz thimbles 


Semiclassical limit 
Let us take h < 1 for computing 


(Oln) = PSU [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around zo's imply that 


"ce. st. (Ole) S » O09). 


RHS is 


[ dze SOM"O(z) ~ S Ole). 
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Complex Langevin and Lefschetz thimbles 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


C, > 0. And, d, = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Consequence 


Naive complex Langevin method cannot explain the Silver Blaze 
phenomenon. 
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Complex Langevin and Lefschetz thimbles 


Proposal for modification 
Assume that 


— Ko, R) 2nh e- S62. 


S" (x. 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P)nsupp(P.) = 0. 


E 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, R) 


2h 
T e SGo)/h 


Sz) Og). 


(0(2)) = DF 


o 
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Complex Langevin and Lefschetz thimbles 


Proposal for modification (conti.) 


If so, defining the phase function 


|S" (Ze)| —iIm S(z¢)/h = 
e) (20)/ Xsupp(Po) (2; Z), 


we can compute 


new ._ (P(n, 24)O(25))n 
OUT — BG. E. — 


This new one is now consistent within the semiclassical analysis. 
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Complex Langevin and Lefschetz thimbles 


Numerical test 1 


S(x) = x*/A — x?/2 + iax (Semiclassical limit is a — oo): 


--- exact 


HEH w/o phase 
s| 1 w/ phase 
20 
[any 
-5 
-10 
0 1 2 3 4 
a 


(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
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CL distribution at —0.5 < 4/U < 1.5 


BU —30 p/U=0 


o 


o 
bo 


0.4 


B[Im(z) — (u + U/2)] 


0 
BRe(z) 
Figure: CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles. 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Dec. 7, 2015 @ Saga 37 / 40 


Complex Langevin and Lefschetz thimbles 


Numerical test 2 
One-site Fermi Hubbard model: (Preliminary) 


2:5 


--- exact 
20l HEH w/o phase 
tQ w/ phase 


-0.5 0.0 0.5 1.0 Ls 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. | hope this can attack the Silver Blaze phenomenon. 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Dec. 7, 2015 @ Saga 38 / 40 


Summary & Conclusion 


Summary and conclusion 
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Summary & Conclusion 


Summary and Conclusion 


e Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 

e One-site Hubbard model is a nice toy model to play with the 
sign problem. 

e Destructive interference among complex phases of Lefschetz 
thimbles play a pivotal role for the baryon Silver Blaze. 

e Naive complex Langevin method is incorrect for generic cases 
and cannot explain baryon Silver Blaze. 

e Modified complex Langevin may become a multi-thimble 
simulator, which has not been achieved in any other algorithms. 
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